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ABSTRACT 

We present the results of a 21-cm HI survey of 27 local massive gas-rich late-stage mergers 
and merger remnants with the Green Bank Telescope (GBT). These remnants were selected from 
the Quasar/ULIRG Evolution Study (QUEST) sample of ultraluminous infrared galaxies (ULIRGs; 
Ls-iooo fim > 10 12 Lq) and quasars; our targets are all bolometrically dominated by active galactic 
nuclei (AGN) and sample the later phases of the proposed ULIRG-to-quasar evolutionary sequence. 

We find the prevalence of HI absorption (emission) to be 100% (29%) in ULIRGs with HI detections, 

100% (88%) in FIR-strong quasars, and 63% (100%) in FIR-weak quasars. The absorption features 
are associated with powerful neutral outflows that change from being mainly driven by star formation 
in ULIRGs to being driven by the AGN in the quasars. These outflows have velocities that exceed 
1500 km s _1 in some cases. Unexpectedly, we find polarization-dependent HI absorption in 57% of our 
spectra (88% and 63% of the FIR-strong and FIR-weak quasars, respectively). We attribute this result 
to absorption of polarized continuum emission from these sources by foreground HI clouds. About 
60% of the quasars displaying polarized spectra are radio-loud, far higher than the ~10% observed 
in the general AGN population. This discrepancy suggests that radio jets play an important role in 
shaping the environments in these galaxies. These systems may represent a transition phase in the 
evolution of gas-rich mergers into “mature” radio galaxies. 

Subject headings: galaxies: active — galaxies: mergers — ISM: jets and outflows 


1. INTRODUCTION 

Recent theoretical work suggests that feedback from 
AGN plays an important role in establishing the present- 
day appearance of massive galaxies. Major galaxy merg- 
ers at high redshifts may trigger powerful starbursts, lead 
to the formation of elliptical galaxies, and account for the 
growth of supermassive black holes (e.g., Sanders et al. 
1988; di Matteo et al. 2005; Springel et al. 2005; Hopkins 
et al. 2005, 2008). In this merger-driven evolutionary 
scenario, two gas-rich progenitors collide to form a com- 
pletely obscured ULIRG where substantial inflows of gas 
trigger starbursts and feed the central black hole(s). As 
the system evolves, the obscuring gas and dust begin to 
disperse, bringing about the formation of a dusty IR- 
excess quasar and eventually a “naked” optical quasar. 
Powerful winds, driven by either the starburst or the cen- 
tral black hole, seem to be required to quench star for- 
mation and to stop the growth of the black hole in order 
to explain the observed Mbh — cr relation (e.g., Murray 
et al. 2005). 

There is growing support for the existence of this type 
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of “negative feedback” in luminous late-stage mergers. 
Most galaxies with star formation rate (SFR) densities 
>0.1 M 0 yr _1 kpc -2 show signs of outflow (e.g., Heck- 
man 2002; Veilleux et al. 2005). These starburst-driven 
winds are often seen in ULIRGs, with mass outflow rates 
of ^0.1-5 times the SFR and velocities ~100-400 km s _1 
(Rupke et al. 2002, 2005a, b,c; Martin 2005, 2006). In 
late-stage mergers and pure Seyfert galaxies, however, 
these winds are sometimes driven by the AGN, reach- 
ing outflow velocities >1000 km s _1 (Rupke et al. 2005c; 
Krug et al. 2010). Similar AGN-driven outflows are also 
found in recently (re)started radio galaxies (e.g., Oost- 
erloo et al. 2000; Morganti et al. 2005b). Several stud- 
ies have also shown that these outflows are often multi- 
phased, with powerful molecular, neutral, and ionized 
outflows having all been detected in the same object (e.g., 
Morganti et al. 2007; Fischer et al. 2010; Feruglio et al. 
2010; Rupke & Veilleux 2011; Sturm et al. 2011). 

In the present paper, we focus on the HI properties of 
the QUEST ULIRGs and quasars, searching for the pres- 
ence and the driving source (AGN or starburst) of neu- 
tral outflows in merging systems as a function of merger 
stage. This paper is organized as follows. In Section 2, 
we discuss how our HI sample is selected. In Section 3, 
we describe our observational setup, data reduction, and 
spectral analysis. We then examine the HI profiles of 
our targets in Section 4 and the detection of polarized 
absorption (including a more sophisticated spectral de- 
composition that boosts our number of absorption de- 
tections) in Section 5. Finally, in Section 6, we discuss 
the possible origins of the HI outflows. We conclude in 
Section 7 with a summary of our findings. Throughout 
this paper, we adopt H 0 = 71 km s _1 Mpc -1 , f l M = 


2 


Teng, Veilleux, & Baker 


0.27, and Da = 0.73 (Hinshaw et al. 2009). 

2. SAMPLE SELECTION 

In order to assess the physical conditions and the dy- 
namical evolution of the neutral gas in late-stage merger 
systems, we selected targets from the QUEST sample of 
ULIRGs and quasars that are bolometrically dominated 
by AGN. The details of the QUEST sample selection are 
outlined in previous QUEST papers (e.g., Schweitzer et 
al. 2006; Netzer et al. 2007; Veilleux et al. 2009a). Briefly, 
the original QUEST sample of ULIRGs begins with the 1- 
Jy sample, a flux-limited sample of 118 ULIRGs selected 
from a redshift survey of the IRAS faint source catalog 
(Kim & Sanders 1998). To these, five well-studied local 
ULIRGs from the Revised Bright Galaxy Survey (RBGS; 
Sanders et al. 2003) were added given the availability of 
archival data from multiple observatories (see Veilleux et 
al. 2009a; Veilleux 2012, for details). The original quasar 
sample contains 25 quasars, including 24 Palomar-Green 
(PG) quasars from the Bright Quasar Sample (Schmidt 
& Green 1983) and B2 2201+31A which satisfies the B 
magnitude criterion of Schmidt & Green (1983). An ad- 
ditional nine PG quasars were added to the sample from 
archival data. The selected quasars fall at the low end of 
the quasar R-band luminosity function. 

The ULIRGs and quasars in the QUEST sample are 
well-matched in redshift (z < 0.3) and in near-IR imag- 
ing show tidal features that indicate merger origins 
(Veilleux et al. 2006, 2009b). The entire QUEST sam- 
ple spans the merger evolution sequence from binary 
starburst-dominated ULIRGs, to completely coalesced 
AGN-dominated ULIRGs, to dust-enshrouded quasars, 
to optical quasars whose dust has been cleared. 

Our HI sample is a subset of the QUEST targets. To 
ensure quality spectra, we require z < 0.12; our 27 tar- 
gets evenly sample coalesced ULIRGs 8 (9 targets), far- 
infr ared (FIR) -strong quasars (/ 6 o/xm//i 5 Mm > 1; 10 tar- 
gets), and FIR- weak quasars (8 targets). The two classes 
of quasars have similar spectral energy distributions rep- 
resentative of AGN, with the FIR-strong quasars having 
an additional contribution to their FIR emission from 
intense star formation (Netzer et al. 2007). The mor- 
phologies of the quasar hosts are in agreement with the 
picture that FIR-weak quasars are at a later evolution- 
ary stage than FIR-strong quasars, following the (near) 
cessation of star formation (Veilleux et al. 2009b). Thus, 
our targets sample the later stages of merger evolution 
from coalescence to the aftermath of significant star for- 
mation. Table 1 summarizes the properties of our HI 
sample. 

3. OBSERVATIONS AND DATA 
3.1. Observational Setup 

The targets were observed with the GBT Spectrometer 
in the L band between 2011 August and 2012 March dur- 
ing the ~160 hours of time allocated for GBT programs 
11B075 and 12A075 (PI: Teng). The data were collected 
in 10 minute on-off source pairs (5 minutes per posi- 
tion) with 12.5 MHz bandwidth, two spectral windows 
centered at the same frequency (1420.4058 MHz), two 

8 This category includes Mrk 273, the only source in our sample 

to have binary nuclei (see Appendix). 


linear polarization modes (XX and YY), and nine sam- 
pling levels. To minimize the effects of radio frequency 
interference (RFI) on the data, we used the notch filter 
for sources with heliocentric recessional velocities below 
17000 km s -1 . The total on-source time for each galaxy 
was determined during the observing runs based on the 
strength of the spectral features relative to the noise or 
up to a maximum of three hours. Most sources were ob- 
served over multiple observing sessions. The circumpolar 
starburst galaxy M82 was observed at the beginning of 
each observing session as a check for instrumental varia- 
tions. 

We also measured the continuum flux densities of our 
sample at 1.42 GHz using the AutoPeak procedure in 
order to check for time variability. However, due to 
the plethora of continuum sources within the large GBT 
beam, the 5 -a confusion limit is ^100 mJy at this fre- 
quency. Therefore, we were only able to obtain confi- 
dent continuum measurements for five targets. These 
values are in agreement with older values obtained with 
the Very Large Array (VLA) for their nuclear core flux 
densities (Condon et al. 2002), despite the much larger 
beam of the GBT. For the remainder of sample, we use 
published values of the core continuum flux densities as 
measured with the VLA (e.g., Barvainis et al. 1996; Con- 
don et al. 1998; Rafter et al. 2009). 

3.2. Data Reduction 

The spectra were reduced using GBTIDL (Marganian 
et al. 2006). Individual five-second records showing large 
broadband harmonic RFI were flagged, and persistent 
RFI spikes that appear in the data were replaced by in- 
terpolation across the affected channels for both spectral 
windows and polarizations. For each polarization, the 
data from the two spectral windows were averaged to- 
gether, and a polynomial of order <5 was fitted over a 
range of ~1000 km s -1 on each side of a galaxy’s systemic 
redshift in the Hanning-smoothed and then decimated 9 
spectrum to subtract the baseline. The same baseline- 
fitting regions and order of polynomial fit were used for 
both polarizations. 

Flux calibration was derived from simultaneous ob- 
servations of an internal noise diode. The baseline- 
subtracted and flux-calibrated scans for each polariza- 
tion were accumulated and averaged for each galaxy. A 
fourth order boxcar smoothing was then applied to the 
averaged data with decimation, resulting in a resolution 
of ^6 km s^ 1 per channel. The two polarizations were 
then averaged together to produce the full intensity spec- 
tra shown in Figure 1. We ensured that our flux calibra- 
tion is correct by comparing the full intensity profile of 
Mrk 231 with that measured by the VLA (Carilli et al. 
1998); both profiles are identical. Discussions of individ- 
ual objects in the Appendix include comparisons with 
previous data, whose generally good agreement with our 
data in terms of velocity profile and flux validate our 
data reduction process. Table 1 also provides a list of 
some observational parameters for our galaxies. 

Nearly all of the galaxies were detected. Two galax- 
ies, PG 0804+761 and F15462-0450, had strong RFI 
spikes overlapping with the expected HI frequencies. At 

9 The total number of channels in the spectrum is decreased, or 
decimated, by a factor equal to the order of smoothing performed. 
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the limit of three hours on source, PG 1126-041 and 
F21219-1757 were not detected, and F04103-2838 and 
F15130-1958 were marginally detected. The detections 
of PG 1426+015 and PG 1617+175 are uncertain be- 
cause their putative emission features have widths sim- 
ilar to those of broadness of the baseline ripples before 
baseline subtraction. 

3.3. Spectral Analysis 

The individual HI profiles in the spectra of our sam- 
ple galaxy are complex combinations of absorption and 
emission. For simplicity, we have measured the emis- 
sion and absorption features separately, but we note that 
this approach may slightly alter the final measured val- 
ues of these quantities. More complex decomposition of 
the emission and absorption is only possible for a few 
sources because of the degeneracy between emission and 
absorption, as discussed in Section 5.2. 

For galaxies with significant emission (peak-to-rms > 
3cr), we measure the central velocity (V©), width (W20, 
W50), and total velocity-integrated flux (/hi) using the 
GBTIDL task gmeasure. In these instances, since there 
are no obvious companions enclosed within the GBT 
beam, the total flux is integrated between velocities at 
which the flux density reaches zero. The central velocity 
is the heliocentric velocity measured at the midpoint of 
the 20% flux level. The width of emission is measured 
at the 20% and 50% levels (W 2 o and W 50 , respectively). 
The measurement errors in these quantities are derived 
using the expressions presented in Wei et al. (2010) fol- 
lowing the methodologies of Schneider et al. (1986, 1990) 
and Fouque et al. (1990). The HI column density is 
calculated using expression 12.17 from Rohlfs & Wilson 
(2000) for optically thin radiation. Assuming the GBT 
antenna temperature to flux density gain conversion of 
~2 K Jy" 1 : 

A/hi = 3.64 x 10 18 /hi cm -2 , (1) 

where /hi is the velocity integrated HI flux in units of 
Jy km s -1 . The neutral gas mass is calculated following 
Haynes & Giovanelli (1984): 

Mhi = 2.35 x 10 5 /hi(^) 2 M q , (2) 

ri 0 

where V cor is the recessional velocity corrected for the 
influence of the Virgo Cluster, the Great Attractor, and 
the Shapley Supercluster as calculated by the NASA Ex- 
tragalactic Database (NED). All of these measured and 
derived quantities are listed in Table 2. 

Only a few absorption features are significantly de- 
tected (by the criterion that peak-to-rms > 3cr). How- 
ever, many of the absorption features are very broad but 
shallow, so the peak-to-rms ratio is not necessarily the 
best determinant of a detection. We also consider a fea- 
ture to be detected if its full width at half maximum 
(FWHM) spans at least 10 velocity channels and it has a 
peak-to-rms > lcr. These absorption features are quan- 
tified using the GBTIDL task fitgauss, with one Gaus- 
sian component at a time added to a spectrum until the 
residuals are statistically featureless. With the excep- 
tion of PG 1119+120, where a narrow absorption feature 
is present within the emission profile, only the absorp- 
tion features are modeled. For PG 1119+120, the emis- 


sion was first modeled with Gaussians and subtracted 
from the data to define a baseline. The absorption was 
then modeled from the emission-subtracted spectrum. 
The Gaussian models are plotted (as components and 
as sums) in Figure 1. The GBTIDL task measures the 
central velocity (V©), FWHM, the “peak” flux density, 
and their measurement errors for each Gaussian compo- 
nent. The full velocity-integrated absorbed flux (/ a bs), 
optical depth (r), and column density (A/n) are then 
calculated from the best-fit models. We assume the HI 
absorbing screen covers the nuclear continuum source, in 
which case the optical depth profile is: 

= ( 3 ) 

*1.4 GHz 

where Shi ( v ) is the HI line profile as a function of velocity 
and Si. 4 ©Hz is the continuum flux density (Whiteoak & 
Gardner 1972). The column density is then calculated 
as: 

^ = 1.823 x 10 18 [ r(u)du cm" 2 , (4) 

S J 

where T s is the spin temperature. In this paper, we will 
assume a spin temperature of 1000 K (see Section 6.1). 
For absorption features in PG 1119+120, PG 1244+026, 
and PG 1440+356, the absorption troughs are deeper 
than their nuclear continua, suggesting that sources of 
continuum emission extend beyond the nucleus. Table 3 
summarizes the measured and derived values for all of 
the detected absorption features. 

As a reality check, we inspected the data to confirm 
that the choice of polynomial order for the baseline fitting 
does not affect the detection of absorption or emission. 
Specifically, we examined whether polynomial orders of 3 
or 5 can introduce ripples into the data that masquerade 
as absorption or emission features. After repeating the 
baseline subtraction for our data using a range of polyno- 
mial orders (1-5), we confirm that the same spectral fea- 
tures are still detected, even for marginal detections such 
as that in F04103-2838. The measured widths and veloc- 
ities are in agreement with those measured with the best- 
fit baseline to within the errors. Therefore, we are confi- 
dent that even the weak detections noted as ’’marginal” 
in Tables 2 and 3 are real. 

4. HI PROFILES OF LATE-STAGE MERGERS 

The individual HI profiles of our targets are discussed 
in the Appendix. In general, we see HI trends along the 
merger sequence but no correlation with infrared lumi- 
nosity. 

All of the detected ULIRGs show HI absorption, typi- 
cally with multiple velocity components, many of which 
are narrow (FWHM < 100 km s -1 ) and some of which 
are redshifted. Only 2/7 (29%) of the detected ULIRGs 
have HI emission. The different velocity absorption com- 
ponents have different origins and are not necessarily 
caused by outflows or infall; kinematics of circumnuclear 
gas disks may mimic the properties of outflows or in- 
fall in single dish data. In a couple of cases, such as 
UGC 05101 and Mrk 273, the absorption profiles are 
nearly symmetrical about zero velocity. For the Type 2 
AGN UGC 05101, the profile is likely from absorption 
due to a slightly inclined rotating HI disk superposed on 
an obscured nucleus. This interpretation is consistent 
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with CO maps of UGC 05101 showing smooth, regular 
rotation of the molecular gas (Genzel et al. 1998; Wil- 
son et al. 2008). On the other hand, the CO velocity 
field of Mrk 273 is very distorted (Downes & Solomon 
1998; Wilson et al. 2008), so it is unlikely that the gas 
is in well-ordered disk-like rotation. The double absorp- 
tion in HI seen in Mrk 273 probably corresponds to the 
binary obscured nuclei in this object. Thus, CO maps 
can help eliminate gas disks as the origin of the absorp- 
tion features. However, only a few galaxies in our sample 
have readily available CO maps so direct comparisons are 
difficult. Comparison of single dish CO and HI profiles 
for the QUEST targets is on-going, but it is beyond the 
scope of the present paper. The more complex profiles, 
such as those seen in F15250+3608, F05189-2524, and 
F15130-1938, may be interpreted as multiple infalling or 
outflowing components. They may also result from the 
superposition of an emission profile and a broad absorp- 
tion profile. The latter case is discussed in more detail 
in Section 5.2. 

In contrast, the HI profiles of the FIR-strong quasars 
are generally detected not only in absorption (88% or 7 of 
the 8 detected sources) but also in emission (75%). The 
absorption features are on average ^60 km s” 1 broader 
than those seen in ULIRGs or approximately three times 
greater than the typical FWHM measurement error in 
the ULIRGs. Both redshifted and blueshifted absorp- 
tion features are seen, although the systemic velocities 
of these systems are uncertain. Optical emission lines 
may underestimate recessional velocity since some frac- 
tion of the line emission may arise from outflowing ma- 
terial. In the case of PG 0050+124 (I Zw 1), this effect 
causes a shift of ^700 km s -1 (see Appendix). A simi- 
lar situation could also hold for PG 1440+356. In three 
cases, PG 1119+120, PG 1244+026, and PG 1440+356, 
the minima of the absorption troughs exceed the nu- 
clear continuum flux densities. This discrepancy implies 
that off-nuclear emission, not accounted for in the nu- 
clear fluxes of Table 1, also contributes to the contin- 
uum emission in these four sources. Given resolution 
limit of the VLA core measurements of PG 1119+120 
and PG 1440+356, the off-nuclear emission must be at 
least ~20-30 kpc from the core; for PG 1244+026, the 
separation from the core must be at least 3 kpc. While 
we cannot formally rule out the possibility that these ab- 
sorption features are absorption of HI against continuum 
from background quasars, this possibility is statistically 
unlikely. Another possibility is that the flux densities of 
the cores in these galaxies are time variable as exempli- 
fied by the variability of the core in PG 1351+640 (see 
Appendix). 

The HI profiles of the FIR-weak quasars are different 
from those of the ULIRGs and FIR-strong quasars. HI is 
seen mostly in emission (88% or 7 of 8 detected sources), 
and only 25% (2/8) show broad absorption. The emission 
profiles are very broad (W$o >500 km s -1 ) and bumpy, 
perhaps beginning to resemble those of regularly rotating 
disks that typically display Gaussian or double-horned 
profiles; the measured widths and velocities are consis- 
tent with previous measurements of quasars (e.g., Bothun 
et al. 1984; Condon et al. 1985; Hutchings et al. 1987; Ho 
et al. 2008; Konig et al. 2009). These properties will be 
revisited in Section 5.2. 


Considering now our entire HI sample, the distribu- 
tion of equivalent widths (EW; defined as the integrated 
flux density for each component divided by the contin- 
uum flux density of the core) is plotted in Figure 2. For 
this definition, the equivalent widths of the absorption 
components are negative. It is clear that the strengths 
of the emission and absorption features progress along 
the merger sequence, with ULIRGs having the weakest 
features, the FIR-strong quasars having the strongest 
absorption features, and the FIR-weak quasars having 
the strongest emission features. In Figure 3, we com- 
pare the FWHMs of emission and absorption with the 
star formation rates of their host galaxies. For emis- 
sion, the FIR-weak sources have the broadest widths, 
the ULIRGs have the narrowest widths, and the FIR- 
strong quasars bridge the gap between those two classes. 
In the ULIRGs, there appears to be a weak correlation 
between the width of the absorption components and the 
SFR, reminiscent of a similar relationship seen in Na I 
(Rupke et al. 2005b). However, given that multiple ab- 
sorption components are associated with a single system 
(and thus a single SFR value), the validity of the connec- 
tion between the SFR and the FWHM of the absorption 
components in ULIRGs is uncertain. On the other hand, 
absorption FWHM is not at all related to SFR in the 
quasars, further suggesting that the AGN play a role in 
these systems. 

5. POLARIZATION-DEPENDENT HI ABSORPTION 

An unexpected outcome of our GBT program is that 
many of the absorption features show a dependence on 
polarization (Figure 4). Indeed, 13 of our 23 GBT- 
detected targets show some profile difference in their two 
polarizations (1 ULIRG, 7 IR-excess quasars, and 5 opti- 
cal quasars). For GBT observations, it is expected that 
there is a <6% difference between polarizations in terms 
of normalization, but the observed spectral shapes should 
be the same in both polarizations (R. Maddalena, private 
comm.). This expectation is confirmed by the fact that 
our M82 data remain in agreement after baseline sub- 
traction. Polarization variation in our QUEST sample 
is not RFI-related, as we have removed time intervals 
where the data are affected by broadband RFI. Since it 
is improbable that the HI screen itself is polarized, the 
underlying continuum emission must be polarized. 

5.1. Time Variability 

The most dramatic examples of polarized spectra are 
seen in PG 1211+143 and PG 1440+356. Not only are 
these spectra polarized, but they are also time variable 
on timescales of ~30 minutes (Figures 5 and 6). In fact, 
of our radio-detected sample galaxies, time variability 
is only seen in sources that have significant absorption 
and spectra that are significantly polarized. A possi- 
ble concern is that since we observed these sources over 
several sessions, variability could be an instrumental ef- 
fect. However, comparing the different sessions’ data for 
M82, we find that time variability and polarization are 
not seen beyond what is expected of typical GBT obser- 
vations. Moreover, session-to-session instrumental vari- 
ability would not explain the 30-minute variability seen 
within the same observing session. Time varying polar- 
izations are not seen only in the targets with the brightest 
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continua, so the variations are not likely to be due to in- 
strumental systematics from poor cancellation of baseline 
ripples in the position-switching scans (O’Neil 2002). 

A different sort of instrumental effect may be leakage 
of linear polarization from the linearly polarized back- 
ground continuum into the total intensity of the absorp- 
tion profiles. We cannot yet exclude this scenario since 
these observations were not designed for polarimetry and 
our current data lack the cross-correlated products (XY 
and YX) that would be needed to confirm these time 
variations. 

If these time variations are real, they could be caused 
by one of several processes. First, they could be due to 
variations in the compact continuum sources behind the 
HI screen. As the X-ray light curves of PG 1211+143 
and PG 1440+356 show (Figure 7), the X-ray emitting 
central region can be highly variable on a timescale of 
about 30 minutes. Although not contemporaneous with 
our GBT data, these X-ray data suggest that the AGN 
continuum can vary on a 30- minute time scale if the radio 
continuum is indeed related to the X-ray producing re- 
gions. However, if the polarization-dependent variations 
are simply a result of the changing compact source con- 
tinuum, then we would expect the continuum to vary in 
block and would observe a change in the equivalent width 
of the HI absorption, not the frequency-dependent vari- 
ations we see in these targets. Second, variability could 
be a result of changes in the structure of the HI screen. 
However, the timescale for the destruction of small HI 
clouds in the interstellar medium is on the order of a 
Myr, depending on the exact properties of the interstel- 
lar medium causing the absorption (e.g., McKee & Cowie 
1977; Nagashima et al. 2006); it is therefore unlikely that 
cloud destruction causes the observed short-term varia- 
tions. Similarly, it is unlikely that the movement of HI 
clouds into and out of our line of sight to the continuum 
source can cause the very broad variations observed in 
these objects. 

The time variability of the polarization-dependent HI 
absorption is easier to explain if the polarized contin- 
uum emission comes from both the nuclear source and an 
extended source (e.g., radio jets). In this picture (Fig- 
ure 8), the HI screen is also extended, covering both the 
nuclear source and extended jets with kinematically dis- 
tinct clouds. As the central source fluctuates, the ratio 
of central-to-jet emission also changes, resulting in the 
observed wavelength-dependent polarization variations. 
In this scenario, there is no need for the structure of the 
HI screen to change on a fast timescale. 

5.2. Decomposing the Spectrum 

Regardless of whether the time variability is real, it is 
certainly the case that the spectra shown in Figure 4 are 
polarized. Again, we checked for baseline-subtraction- 
related effects on the detection of polarized absorption 
and found that the choice of the polynomial order of 
the baseline fit matters little to the detection of these 
polarization-dependent absorption features. The most 
obvious example of polarization not due to observational 
systematics is seen in the spectra of PG 0050+124 (I Zw 
1). This HI profile is not polarized at relative velocities 
of 500-1000 km s -1 but becomes polarization-dependent 
between relative velocities -500 and +500 km s -1 (Fig- 
ure 4). This behavior argues against attribution of ap- 


parent polarization dependencies to poor baseline sub- 
traction or session-to-session variations in our GBT ob- 
servations. At least in this one case, the polarization de- 
pendence must be due to something physically happen- 
ing at relative velocities between -500 and +500 km s -1 . 
When we compare the HI profiles of each galaxy in the 
two separate polarizations, it appears that one polariza- 
tion typically has higher flux than the other, with the 
exception of PG 1440+356, in which there is a veloc- 
ity shift instead. The effect is not systematic in the 
sense that one polarization is always higher than the 
other; each is nearly equally likely to be the polarization 
with the higher intensity. In some cases, particularly in 
PG 1351+640 and PG 1411+442, the lower-intensity pro- 
file appears to be the same as the higher-intensity profile 
with the addition of a broad absorption component. 

To follow up on this intriguing hint, we produced the 
difference spectra (the spectrum in the polarization with 
the more negative HI profile minus the spectrum in the 
polarization with the more positive HI profile) shown 
in Figure 9. The absorption features in these spectra 
are then modeled the same way as before, using the 
GBTIDL task fitgauss. When these absorption mod- 
els are included, many of the full intensity spectra that 
need multiple components to explain the blue- and red- 
shifted components can be explained simply by HI emis- 
sion that is combined with a single broad absorption com- 
ponent. This pattern is most obvious in PG 1351+640 
and PG 1411+442. The results of this type of complex 
decomposition of the spectra are provided in Tables 4 and 
5. Similar decompositions can explain the shapes of the 
full intensity spectra of the ULIRGs F05189-2524 and 
F15130-1958, and the FIR-strong quasar PG 1244+026. 
However, since those spectra are not polarized, we can- 
not create difference spectra to model the broad absorp- 
tion component. The fitgauss task statistically prefers 
the multiple-component models to the single-absorption 
model for these spectra. Overall, with the exception of 
PG 1501+106, the decompositions result in stronger fea- 
tures in both absorption and emission. 

Given the results of this exercise, we find the instances 
of HI absorption (emission) to be 100% (29%) in detected 
ULIRGs, 100% (88%) in FIR-strong quasars, and 63% 
(100%) in FIR-weak quasars. The more complex de- 
composition of the spectra does not drastically change 
the distribution of EWs (Figure 10) compared to Fig- 
ure 2. Figure 11 provides a comparison of the FWHMs 
of the emission and absorption features from the complex 
modeling with the star formation rates of these galaxies. 
As in Figure 3, the FIR-weak quasars have the broad- 
est emission features and the ULIRGs have the narrow- 
est emission features. In terms of absorption, there is 
no obvious relationship between the absorption feature 
widths and the star formation rates. As before, the FIR- 
weak quasars have the broadest absorption features, with 
the FIR-strong quasars spanning the range of widths be- 
tween FIR-weak quasars and ULIRGs. 

6. THE PROPERTIES OF HI OUTFLOWS 

As discussed in Section 4, the absorption components 
are not necessarily caused by infall or outflow. It is pos- 
sible, as seen in UGC 05101 (edge on) and Mrk 231 (face 
on), that these features can be associated with the kine- 
matics of gas disks. A circumnuclear gas disk seen nearly 
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edge on with a non-isotropic gas mass (or background 
continuum) distribution can mimic characteristics of in- 
fall and outflow. Unfortunately, our single dish data can- 
not unambiguously exclude gas disks as the source of the 
features we observe, but it is unlikely that all of our tar- 
gets with absorption features have gas disks that are seen 
precisely edge-on. 

Keeping this caveat in mind, we assume in the follow- 
ing discussion an outflow origin only for blueshifted HI 
absorption features with Hi-corrected heliocentric veloc- 
ities (V 0i hi 10 ) below -50 km s -1 (Rupke et al. 2005a). 
A number of mechanisms have been proposed to explain 
such HI outflows. They include adiabatically expand- 
ing broad emission line clouds (e.g., Elvis et al. 2002), 
starburst- and AGN-driven galactic winds (e.g., Heck- 
man et al. 1990; Veilleux et al. 2005), and AGN-driven 
jets interacting with the surrounding interstellar medium 
(ISM; e.g., Oosterloo et al. 2000). These mechanisms pre- 
dict different outflow locations (from parsec to kiloparsec 
scales) and geometries (wide angled vs. collimated). 

6.1. Mass Outflow Rates 

Following Heckman et al. (2000) and Rupke et al. 
(2002), the mass outflowing (M) at a velocity v is calcu- 
lated using: 

— = 20( — )( r * )( ~^ HI )( — ) 

Mq yr- 1 M-7T 1 kpc'nO 21 cm -2 ^200 km s -1 ’ 

( 5 ) 

assuming spherical geometry, mass conservation, and a 
wind that flows from a distance of r* to infinity. The out- 
flow is also assumed to subtend a solid angle of H within 
which the gas has a covering fraction of 1. To facilitate 
comparisons with the values derived from young radio 
galaxies by Morganti et al. (2005b), we use the same as- 
sumptions of If = 7r and r* = 1 kpc. Applying these as- 
sumptions and Equation 5 to the blueshifted absorption 
components where v = Vq, hi, we find mass outflow rates 
M (wide) up to more than 6000 Mq yr -1 (summarized 
in Table 6). These values are much larger than those of 
<60 Mq yr -1 measured by Morganti et al. (2005b) and 
the range measured in ULIRGs (13-133 M 0 yr -1 ; Rupke 
et al. 2002). 

Many of the mass outflow rates we infer seem unrealis- 
tically large; such high rates have not been measured in 
other astrophysical objects. These values do, however, 
depend on the assumptions we have made above. The 
mass outflow rate depends linearly on the size scale of 
the outflow region and the derived HI column density. 
Recall from Equation 4 that the column density depends 
on the spin temperature and optical depth profile, which 
in turn depends on the brightness of the background con- 
tinuum. There are thus three key factors that affect our 
estimates of a mass outflow rate: the spin temperature 

10 These systemic velocities are relative to the HI emission rather 
than the optical emission. See Table 6 for more details. A remain- 
ing source of uncertainty in the interpretation of HI line profiles is 
the question of whether the midpoint of the 20% flux level is always 
a good approximation of the true systemic velocity. Emission from 
molecular lines is a more reliable tracer of systemic velocity than 
HI, so the ideal solution here would be acquisition of single-dish 
CO spectroscopy for the full QUEST sample. Unfortunately these 
observations are still pending. 


of the gas, the assumed distance, and the background 
continuum. The assumption of 1000 K for the HI spin 
temperature seems to be reasonably conservative. In the 
vicinity of powerful AGN, the extreme conditions of the 
gas and the presence of shocks can result in spin temper- 
atures >3000 K (e.g., Holt et al. 2006), although we note 
that absorption-line opacity will be dominated by the 
coldest material in any multiphase medium. Spatially 
resolved interferometric maps of several quasars suggest 
that the outflow regions have size scales ^0.2-1. 6 kpc 
(Oosterloo et al. 2000; Morganti et al. 2003; Emonts et 
al. 2005; Morganti et al. 2005a). Adopting the smaller 
distance of 0.2 kpc instead of 1 kpc would reduce the de- 
rived mass outflow rates by a factor of 5, more consistent 
with those derived for ULIRGs. 

Finally, the estimated HI column density also depends 
strongly on the geometry of the source of continuum 
emission. We have assumed that this emission comes 
entirely from the core. However, the absorption troughs 
in several quasars in our sample are deeper than the core 
fluxes. This excess implies that off- nuclear emission also 
contributes to the observed continuum. In this case, r in 
Equation (3) is overestimated and so is A r pn , by a factor 
equal to the actual nuclear/total continuum flux ratio, 
although a larger source size may compensate to some 
extent. 

6.2. Jet-induced Outflows 

If the off-nuclear continuum source is indeed physically 
related to the target and not a background quasar that 
happens to lie within the GBT beam, our polarized spec- 
tra give more credence to the idea that radio jets play 
an important role in driving HI gas outflows (e.g., Mor- 
ganti 2011). Several sources (see Appendix) in our sam- 
ple with polarized spectra have spatially resolved radio 
jets detected at higher frequencies. The resolutions of 
these VLA data imply that the jets are launched within 
~0.4-0.9 kpc of the core. Considering a jet-induced out- 
flow scenario for the HI absorption features, we again 
apply Equation 5 to the blueshifted outflow components 
and assume r* = 0.5 kpc and a narrow opening angle of 
5° for the jet. Under these new assumptions, the range 
of mass outflow rates is 0.2-31 M 0 yr -1 (see col. (8) 
of Table 6). We consider these values lower bounds on 
the actual mass outflow rates, while those derived using 
H = 7r and r» = 1 kpc (col. (7) of Table 6) are upper 
bounds. 

Under the jet scenario, the kinetic energy release rate 
can be calculated using the equation from Heckman et 
al. (2000): 

^ — 2.7x 10 41 ( — If r * )( Nhi V V 

erg s — 1 4 tt 1 kpc 10 21 cm -2 zOO km s -1 

( 6 ) 

The kinetic energy in the outflow is on the order of 10 38 
to 10 42 erg s _1 for our sample. 

The association of jets with the detected HI absorption 
features is also seen in the distribution of radio loud- 
ness of our sample. Only one quasar in our sample is 
considered radio-loud based on optical-to-radio flux ra- 
tios (see Appendix), but optical data may be subject to 
contamination from the host galaxy, resulting in an un- 
derestimate of radio loudness. Using the 1.4 GHz to 2- 
10 keV ratio as an indicator of radio loudness and a radio- 
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loud cutoff value empirically derived from 1600 AGN (La 
Franca et al. 2010), we find that many objects in our 
sample are radio-loud (Table 1). These results are con- 
sistent with those found in a recent study by Chandola et 
al. (2012) that suggest higher detection rates of HI ab- 
sorption toward extended compact steep-spectrum and 
gigahertz-peaked spectrum sources than toward the cores 
of radio galaxies. The environments in which the jets in- 
teract with the interstellar medium can be complex, and 
we may be seeing a strongly polarized continuum being 
absorbed selectively by multiple neutral hydrogen clouds, 
resulting in polarized HI absorption. 

The fraction of radio-loud sources with associated HI 
absorption (58%) is far higher than the ~10% observed 
in the general population of AGN. Given the high star 
formation rates of the ULIRGs, the circumnuclear star- 
burst can contaminate the 1.4 GHz measurements. Using 
the 1.4 GHz to SFR correlation of Bell (2003), we cor- 
rected the 1.4 GHz flux for star formation. Of the three 
ULIRGs with blueshifted absorption and X-ray detec- 
tions, all (100%) are radio-loud. Similarly, of the six 
FIR-strong quasars with detected HI outflows and X-ray 
emission, three (50%) are considered radio-loud; one of 
the two (50%) FIR-weak quasars with HI outflows and X- 
ray flux are radio-loud. This pattern recalls an idea put 
forward by Wilson & Colbert (1995) on the formation of 
radio galaxies to explain why radio-quiet AGN outnum- 
ber radio-loud AGN. The authors suggested that major 
mergers, such as the progenitors of our sample of galax- 
ies, result in coalesced supermassive black holes with high 
spin values due to the conservation of angular momen- 
tum. Substantial energy is then extracted from these 
fast-spinning black holes to drive jets, yielding radio- 
loud AGN. Thus, younger radio galaxies may have more 
powerful jets and therefore stronger jet-induced outflows. 
This scenario is consistent with the observations of Mor- 
ganti et al. (2005b), whose sample of young or recently 
restarted radio galaxies show similar HI outflows. The 
Wilson & Colbert (1995) picture is also consistent with 
theoretical models of Tchekhovskoy et al. (2010) that 
suggest the radio loud/quiet dichotomy in AGN is due 
to the spin of the central black hole. 

6.3. HI Outflows Along the Merger Sequence 

Even though our sample with mass outflow rates is 
relatively small, we see some trends in our outflow mea- 
surements (Table 6 and Figure 12). However, given the 
small number of detected HI outflows, particularly with 
only two measurements for the FIR-weak quasars, we 
cannot statically confirm that the distributions we see in 
Figure 12 is significant. There is also no obvious rela- 
tionships between mass outflow rate and the bolometric 
luminosity, AGN luminosity, or radio loudness for our 
sample. 

Systems with deep absorption and outflows become 
more prevalent as we progress along the merger sequence 
from the newly coalesced ULIRGs to the FIR-strong 
quasars. The trend between outflow kinematics and star 
formation rates among the ULIRGs (Figures 3 and 11) 
suggests that the HI outflows in these objects are driven 
mainly by the starbursts, with a median mass outflow 
rate in the jet scenario (M (jet) ~1 Mq yr -1 ) similar 
to their SFRs and previous Na I-based outflow measure- 


ments (Rupke et al. 2002, 2005a, b,c; Martin 2005, 2006). 
ULIRGs also tend to have the narrowest absorption fea- 
tures. In the FIR-strong quasar systems, the median 
mass outflow rate is M (jet) ~5.8 Mq yr -1 , considerably 
larger than their median SFR of ~1.0 Mq yr -1 . While 
the median mass outflow rate of the FIR-weak quasars 
is lower, M (jet) ~0.3 Mq yr -1 , so is their median SFR 
(~0.6 Mq yr -1 ). Thus, the mass-loading factor, i.e., 
the ratio of mass outflow rate to SFR, is much greater 
in the FIR-strong than the value of near unity derived 
in ULIRGs. The same value is much lower in the FIR- 
weak quasars, but this is difficult to compare given that 
this is based on a single measurement. However, the 
distinction between the ULIRGs and FIR-strong quasars 
further supports the idea that the AGN is playing a role 
in driving HI outflows in the quasars before the end of 
large-scale star formation. 

HI in emission is much more prevalent among the 
quasars than the ULIRGs, and the HI emission profiles 
increasingly resemble those of regularly-rotating neu- 
tral gas disks as we move from FIR-strong to FIR- 
weak quasars. The amount of HI seen in emission in- 
creases modestly from ~4xl0 10 Mq among the FIR- 
strong quasars to ~7xl0 10 Mq among the FIR-weak 
quasars. These masses are consistent with those mea- 
sured in low-redshift quasars (Konig et al. 2009) and the 
high end of the gas mass range seen in nearby active 
galaxies (Ho et al. 2008). We may be witnessing the cre- 
ation of HI disks from residual neutral gas at the centers 
of merger remnants. 

7. SUMMARY 

We have analyzed of GBT HI spectra of 27 ULIRGs 
and quasars selected from the QUEST sample. We find 
that: 

1. As mergers evolve from fully coalesced ULIRGs 
to FIR-strong quasars to FIR-weak quasars, the 
characteristics of their HI profiles also change from 
mostly absorption in the ULIRGs to mostly emis- 
sion in the FIR-weak quasars, with the FIR-bright 
quasars having intermediate properties. The kine- 
matics of these systems change from narrow ab- 
sorption features likely due to outflows driven by 
starbursts in ULIRGs, to broad absorption features 
from jet-driven outflows driven by AGN in FIR- 
strong quasars, and finally to the emission profiles 
more characteristic of gas in rotation in FIR-weak 
quasars. 

2. A large fraction of the quasars exhibit polarization- 
dependent absorption features, with 78% of FIR- 
strong quasars and 63% of FIR-weak quasars show- 
ing these features compared with 13% in ULIRGs. 
Some of these features vary on short time scales 
(down to ~30 minutes). Neither time variabil- 
ity nor polarization-dependent spectral features are 
seen in contemporaneous observations of M82 (as 
expected), reducing the likelihood that these ef- 
fects are instrumental artifacts. These features are 
likely due to HI absorption of polarized nuclear and 
off-nuclear continuum emission. The large frac- 
tion (~60%) of radio-loud sources with HI outflow- 
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associated absorption suggests a close relationship 
between radio jets and the HI screen. 


3. The polarization-dependent absorption features are 
most likely due to the emergence of radio jets 
in these systems, contributing to the continuum 
emission against which the absorption features are 
detected. The ~60% radio-loud fraction of our 
sources with HI outflows is far higher than the 
~10% seen in the general population of AGN. This 
distinction suggests that radio jets play an impor- 
tant role in shaping the environments in these sys- 
tems, which may tract a transition phase in the 
evolution of gas-rich mergers into “mature” radio 
galaxies. 
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APPENDIX 

NOTES ON INDIVIDUAL GALAXIES 

In this appendix, we describe the HI spectral properties of our GBT sample. When appropriate, we compare the 

present data with previous observations of these galaxies. 

PG 0007+106: The only radio-loud object in our sample based on classic radio-to-optical criteria. This source is 
classified as a triple system in the third Zwicky catalog (III Zw 2), and the observed velocity range is consistent 
with the A ( v sys = 26921 km s -1 ) and C (v sys = 26981 km s -1 ) components of the system. The GBT detection 
is assumed to be that of the brightest component, A. The HI emission profile of the full intensity spectrum agrees 
with the old Arecibo spectrum of Hutchings et al. (1987). The absorption feature, also noted by Hutchings et 
al. (1987), is more prominent in the XX spectrum of the GBT data. 

PG 0050+124 (I Zw 1): The heliocentric radial velocity listed in Table 1 for this source is from optical spectroscopy 
by Ho & Kim (2009). The optically-derived velocity is ~700 km s -1 lower than the Hl-derived value of 18285 
km s -1 (Bothun et al. 1984). The optical emission lines may arise in part from outflowing gas, so the optically- 
derived value may is underestimate the actual systemic velocity of this object. The shape of the full intensity 
profile (Figure 1) at 18285 km s -1 agrees with those previously published by Bothun et al. (1984) and Hutchings 
et al. (1987), and matches expectations for a regularly rotating disk. However, at lower velocities, the profile 
becomes polarized and a broad absorption feature is centered at the optically-derived velocity (see Figure 9). 

F04103-2838: There is a marginal detection of an emission feature and a weak but broad blueshifted absorption 
feature present. 

F05189-2524: No emission is seen in the data. However, several absorption components are detected. An alternative 
explanation for the spectral shape is the superposition of an emission profile with a broad absorption component. 

PG 0804+761: Strong RFI is present at the expected frequency of this galaxy. It is therefore excluded from our 
analysis. 

PG 0844+249: The broad emission profile is asymmetric. Although the profiles in the two different polarizations 
differ slightly, they have shapes similar to that observed by Ho et al. (2008). Our integrated GBT flux density 
is larger than that measured by Ho et al. (2008). The emission feature redward of the peak may be due to a 
companion at v sys = 19601 km s _1 . 

UGC 05101: Multiple absorption components are seen in the GBT data. We interpret the HI absorption profile as 
a highly absorbed core coincident with an inclined regularly rotating disk. 

PG 1119+120: This object was previously observed by Bothun et al. (1984) and Hutchings et al. (1987) with Arecibo. 
The profiles of the older spectra resemble our GBT observation, although the narrow absorption feature is not 
apparent in the lower-resolution Arecibo data. There is no obvious detection of a companion in the GBT data 
near 16000 km s -1 as suggested by Hutchings et al. (1987). High resolution VLA maps at 4.8 GHz show a bent 
jet extending to the northeast (Leipski et al. 2006). At the resolution of these VLA maps (~0.75"), the jet is 
resolved at a distance of ~0.72 kpc. 

PG 1126-041: This target was statistically undetected after three hours on-source, although there is a very weak 
broad absorption feature centered near its systemic velocity. 

PG 1211+143: This broad absorption line system was undetected by Condon et al. (1985); when it was observed by 
Ho et al. (2008) with Arecibo, standing waves were present in the data due to the strong continuum. 

PG 1229+204: This object was observed by Ho et al. (2008). The profile of the GBT spectrum roughly agrees with 
the earlier data. 

PG 1244+026: A broad, blueshifted absorption feature is present. As for F05189-2524, an alternative explanation 
for the spectral shape may be emission absorbed by a single broad absorption feature. 
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Mrk 231: The deep absorption trough is associated with the face-on disk of the galaxy (Carilli et al. 1998; Ulvestad 
et al. 1999). Our GBT full intensity HI profile of this sources mostly agrees with these previous observations, but 
we detect an additional absorption feature centered at relative velocity ~-300 km s -1 that was first suspected 
by Morganti (2011). This feature is slightly polarized, with the XX spectrum having lower intensity. Although 
the polarization is not as high as those seen in Figure 4, if we give the spectra of Mrk 231 the same treatment 
detailed in Section 5.2, we find that this feature has FWHM ^1000 km s -1 corresponding to a mass outflow rate 
of ~0.3 M 0 yr -1 in the jet scenario. The polarized feature is thus likely associated with the radio jet (Carilli et 
al. 1998; Ulvestad et al. 1999). The velocity of this jet-enhanced outflow is similar to the value measured using 
Na I in the jet-disturbed region (-1400 km s -1 ; Rupke & Veilleux 2011). 

Mrk 273: This source is the only binary in our sample; its nuclear separation is ~750 kpc (Scoville et al. 2000). X-ray 
observations by Iwasawa et al. (2011) identified its southwestern nucleus as the Seyfert 2 nucleus and suggested 
that the northern nucleus may host a heavily obscured AGN. The double-peaked HI absorption profile is very 
symmetric and is likely due to the two obscured nuclei. There appears to be a weak absorption feature in the 
YY spectrum at -500 km s -1 . 

PG 1351+640: The XX spectrum appears to show emission, while YY shows broad absorption. The HI absorption 
data may also be variable. Leipski et al. (2006) note a jet-like structure in their 4.8 GHz VLA maps and time 
variability in the core flux density. The resolution of the Leipski et al. (2006) data (~0.56") corresponds to a 
distance of ~0.91 kpc. Although on a different spatial scale, this structure is consistent (similar position angle) 
with a single-sided extension hinted at by the 8.4 GHz VLB A snapshot data presented by Blundell & Beasley 
(1998). This feature is not seen in deeper VLBA imaging at multiple frequencies, but it should be noted that 
the data from Ulvestad et al. (2005) are only in the left-circular polarization. 

PG 1411+442: The spectrum of this object is very similar to that of PG 1351+640. Emission is seen in the XX 
spectrum and a two-component absorption feature in the YY spectrum. 

PG 1426+015: Ho et al. (2008) detect a fairly narrow (U+o = 357.4kms _1 ) and well-defined HI line centered at 
v sys = 25975.2 km s -1 . We have examined our spectrum at the corresponding frequency and see hints of emission 
over a wider range of velocities; however, this feature has such low significance and is so sensitive to the details 
of baseline fitting we classify its detection status as ” uncertain.” 

PG 1440+356: The profile of the HI absorption feature depends slightly on the polarization. Exceptionally, the 
feature is redshifted by ~300 km s -1 with respect to systemic velocity rather than being blueshifted. A possible 
explanation is that the systemic velocity taken from NED and based on optical emission lines underestimates the 
actual systemic velocity because the line-emitting gas arises in part from outflowing material. Another possibility 
is that we may be detecting infalling HI. Emission is also present in the blue half of the spectrum, although this 
is only seen at one of several epochs. 

PG 1448+273: The weakly polarized spectra show blueshifted absorption. 

PG 1501+106: A weak emission feature is present in the data. A blueshifted absorption feature is also visible, 
particularly in the XX spectrum. A lot of the data for this source were affected by small ripples in the baseline. 
As a reality check, we compared our data with similar GBT data from program 10A-070 (PI: Kannappan). The 
ripples also appear to be present in a large fraction of their data. The persistent ripples could be due to harmonic 
waves from a strong RFI signal at a frequency that is slightly outside the detector frequency range. Due to the 
expected weakness of the features, we conservatively flagged all of the data that showed these ripples. 

F15130-1958: Multiple weak absorption components are present in this object. 

F15250+3608: This is the only ULIRG in our sample that shows distinct polarization-dependent HI absorption. The 
spectra are very similar to those of PG 1351+640 and PG 1411+442, where weak emission is seen in the XX 
spectrum and absorption in YY. 

F15462-0450: Strong RFI is present at the expected observed frequency of this galaxy. It is therefore excluded from 
our analysis. 

PG 1617+175: Broad emission is detected in this object. The profiles of the emission differ slightly in the two 
polarizations. Notably, the emission feature in the YY spectrum seems to be affected by an additional absorption 
feature (Figure 9). 

F21219-1757: This object was not detected. 

PG 2130+099: This target was observed by Ho et al. (2008). Our full intensity profile of this source agrees with the 
older data from Arecibo. A blueshifted absorption feature is also seen in the XX spectrum. Leipski et al. (2006) 
detected what looks like jet-associated hot spots emanating from the core in 4.8 GHz A-array VLA data. The 
resolution of these VLA data (^0.35") implies that the jet can be resolved at a distance of ~0.42 kpc. 

PG 2214+139: The HI spectrum of this object is dominated by a broad emission feature centered at systemic velocity. 
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TABLE 2 

Measurements of HI in Emission 


Source 

(i) 

y © 

(km s 1 ) 
(2) 

W 20 

(km s^ 1 ) 

(3) 

W 50 

(km s x ) 

(4) 

Peak 

(mjy) 

(5) 

S/N 

(6) 

/hi 

(Jy km s — 1 ) 
(7) 

N m 

(10 18 cm -2 ) 
(8) 

M H i 

(io 10 M e ) 
(9) 

PG 0007+106 

-216+ 16 

1063 ± 49 

894 ± 33 

4.30 

5.58 

2.57 + 0.14 

9.35 + 0.49 

8.46 ± 0.45 

PG 0050+123 

473 ± 17 

789 ± 51 

435 ± 34 

7.00 

7.45 

2.43 + 0.14 

8.85 + 0.51 

3.44 + 0.20 

F04103-2838 (M) 

129 ± 21 

151 + 63 

21 ± 42 

3.44 

3.91 

0.17 + 0.06 

0.62 + 0.22 

0.97 + 0.34 

PG 0844+249 

243 + 8 

541 ± 23 

436 ± 15 

8.83 

9.20 

2.32 + 0.12 

8.44 ± 0.43 

4.05 + 0.21 

PG 1119+120 

-18 + 2 

226 + 7 

221 ± 5 

5.35 

6.45 

0.58 ± 0.07 

2.11 + 0.24 

0.65 ± 0.07 

PG 1211+143 

522 + 8 

200 ± 25 

174 + 17 

3.85 

4.18 

0.04 ± 0.07 

0.14 + 0.25 

0.11 + 0.20 

PG 1229+204 

182 ± 11 

600 ± 34 

504 ± 23 

5.37 

5.84 

1.84 + 0.12 

6.70 + 0.43 

3.24 + 0.21 

Mrk 273 

407 ±11 

116 + 32 

62 + 21 

5.94 

4.53 

0.37 + 0.07 

1.35 + 0.26 

0.24 + 0.05 

PG 1351+640 

-205 ± 21 

585 ± 64 

480 ± 43 

2.64 

3.34 

0.82 ± 0.10 

2.98 + 0.38 

2.73 ± 0.34 

PG 1426+015 (U) 

107 + 35 

1131 ± 104 

591 ± 69 

5.12 

4.69 

2.08 ± 0.20 

7.57 + 0.72 

6.86 + 0.65 

PG 1440+356 

-294 + 13 

477 ± 39 

395 ± 26 

3.01 

4.81 

0.48 ± 0.07 

1.73 + 0.27 

1.30 + 0.20 

PG 1448+273 

-50 + 18 

97 + 54 

23 + 36 

2.04 

3.26 

0.11 + 0.03 

0.41 + 0.12 

0.21 ± 0.06 

PG 1501+106 

61 ± 21 

137 ± 63 

22 ± 42 

2.58 

3.37 

0.32 + 0.05 

1.16 + 0.17 

0.20 + 0.03 

PG 1617+175 (U) 

-118+ 15 

1139 + 44 

1057 ± 30 

6.94 

4.36 

3.88 + 0.30 

14.12 ± 1.08 

21.21 ± 1.62 

PG 2130+099 

224 + 9 

599 ± 26 

518 ± 18 

7.86 

6.95 

2.30 + 0.15 

8.37 + 0.53 

3.84 ± 0.24 

PG 2214+139 

104 ± 14 

876 ± 42 

567 ± 28 

11.22 

8.53 

4.43 + 0.21 

16.13 + 0.75 

7.99 ± 0.37 


Note. — Col.(l): Galaxy name where (M) and (U) indicate marginal and uncertain detections, respectively, as noted in 
Section 3.2. Col. (2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table 1, and its 
associated measurement error. Col. (3): Velocity width of the emission measured at the 20% level and its error. Col. (4): Velocity 
width of the emission measured at the 50% level and its error. Col. (5): Maximum intensity of the spectrum. Col. (6): Signal- 
to-noise ratio given by the peak intensity relative to the channel rms listed in Table 1. Col. (7): Velocity-integrated HI flux in 
emission and its error. Col. (8): Derived HI column density calculated using Equation 1. Col. (9): Derived HI gas mass from 
emission and its error calculated using Equation 2. 
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TABLE 3 

Measurements of HI in Absorption 


Source 


FWHM 

“Peak” 

S/N 

N m 

Tmax 

fahs 


(km s *) 

(km s -1 ) 

(mJy) 

(10 : 

22 cm ~ 2 ) 


(Jy km s -1 ) 

(i) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

(6) 

( 7 ) 

(8) 



Strong Detections: peak-to 

-rms > 3(J 




F05189-2524 

160 + 3 

79 + 8 

-2.74 + 0.23 

4.73 

1.50 

0.099 

-0.04 

UGC 05101 

-251 ± 11 

419 ± 26 

-6.69 + 0.22 

4.59 

3.24 

0.040 

-0.54 


8 + 3 

96 + 7 

-9.76 + 0.59 

6.70 

1.09 

0.059 

-0.18 


203 + 9 

219 + 22 

-5.23 + 0.34 

3.59 

1.32 

0.031 

-0.22 

PG 1119+120 

-88 + 3 

22 + 6 

-4.20 + 0.87 

5.06 



-0.02 

PG 1211+143 

-227 ± 1 

460 + 3 

-37.67 + 0.16 

40.95 

5.73 

0.065 

-3.07 

PG 1244+026 

-381 ± 8 

226 ± 20 

-3.07 + 0.23 

3.05 



-0.13 

Mrk 231 

-6 + 2 

199 + 6 

-16.87 + 0.41 

11.07 

2.15 

0.056 

-0.64 

Mrk 273 

-122 + 8 

334 ± 21 

-11.41 + 0.41 

8.71 

5.60 

0.088 

-0.74 


187 + 6 

171 ± 13 

-10.62 + 0.64 

8.11 

2.66 

0.081 

-0.35 

PG 1440+356 

325 + 4 

333 + 8 

-9.57 + 0.21 

15.27 



-0.57 



Weak Detections: peak-to- 

rms < 3cr 




PG 0050+124 

15 ± 10 

53 + 24 

-1.34 + 0.51 

1.42 

1.62 

0.161 

-0.01 

F04103-2838 (M) 

-409 ± 29 

449 ± 70 

-1.17 + 0.16 

1.33 

8.10 

0.094 

-0.09 

F05189-2524 

-374 ± 9 

118 + 22 

-1.29 + 0.20 

2.22 

1.03 

0.046 

-0.03 


-79 ± 13 

163 ± 30 

-1.09 + 0.17 

1.88 

1.20 

0.038 

-0.03 


347 + 11 

127 + 27 

-1.05 + 0.19 

1.81 

0.90 

0.037 

-0.03 

PG 1211+143 

-644 ± 25 

214 ± 60 

-1.14 + 0.25 

1.24 

0.08 

0.002 

-0.04 


766 ± 11 

116 + 26 

-2.47 + 0.48 

2.68 

0.09 

0.004 

-0.05 

PG 1244+026 

183 ± 18 

268 ± 43 

-1.52 + 0.21 

1.52 

59.07 

1.407 

-0.08 

Mrk 231 

-280 ± 9 

96 + 21 

-3.14 + 0.58 

2.06 

0.19 

0.010 

-0.06 

PG 1351+640 

-224 ± 10 

115 ± 25 

-1.53 + 0.28 

1.94 

1.02 

0.046 

-0.03 


336 ± 37 

516 ± 98 

-0.97 + 0.14 

1.23 

2.88 

0.029 

-0.09 

PG 1411+442 

-446 ± 10 

319 + 23 

-2.60 + 0.16 

2.81 



-0.14 


-41 ± 10 

214 + 25 

-1.70 + 0.17 

1.84 

>39.78 

1.139 

-0.06 


319 ± 10 

143 ± 24 

-1.76 + 0.26 

1.90 

>28.08 

1.217 

-0.04 

PG 1448+273 

-441 ± 14 

392 ± 34 

-1.69 + 0.13 

2.70 

55.58 

0.829 

-0.12 


144 ± 19 

195 ± 45 

-0.74 + 0.15 

1.18 

10.27 

0.283 

-0.03 

PG 1501+106 

-229 ± 12 

304 ± 28 

-2.12 + 0.17 

2.77 

>83.12 

1.884 

-0.12 


340 ± 20 

195 ± 48 

-0.98 + 0.21 

1.28 

>17.46 

0.498 

-0.04 

F15130-1958 (M) 

-847 ± 10 

89 + 23 

-2.07 + 0.47 

1.66 

3.86 

0.232 

-0.03 

-451 ± 15 

168 ± 36 

-1.75 + 0.32 

1.40 

6.10 

0.192 

-0.05 


199 ± 16 

182 ± 38 

-1.74 + 0.32 

1.40 

6.56 

0.191 

-0.05 

F15250+3608 

-42 + 8 

83 + 18 

-2.57 + 0.49 

2.21 

2.42 

0.154 

-0.04 


229 + 8 

107+19 

-2.79 + 0.43 

2.40 

3.41 

0.168 

-0.06 


Note. — Col.(l): Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic 
velocity presented in Table 1, as modeled using a Gaussian model, and its associated measurement error. Col. (3): 
Full width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. Col. (4): 
Peak absorption flux and its error as measured using a Gaussian model. Col. (5): Signal-to-noise ratio measured 
as a ratio between the peak flux and channel rms listed in Table 1. Col. (6): Integrated HI column density based 
on the absorption model, assuming a spin temperature of 1000 K and using Equation 4 based on the optical depth 
profile derived using Equation 3. For PG 1411+442 and PG 1501+106, the 1.4 GHz core flux densities are only 
upper limits, so their derived JVhi values are lower limits. Col. (7): Peak optical depth based on Equation 3. Col. (8): 
Velocity-integrated HI flux in absorption measured from the best-fit Gaussian component models. 
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TABLE 4 

Measurements of HI in Emission from Complex Spectral Modeling 


Source 

(i) 

v 0 

(km s - 1 ) 
(2) 

W 20 

(km s -1 ) 

(3) 

W 50 

(km s 1 ) 

(4) 

“Peak” 

(mJy) 

(5) 

S/N 

(6) 

/hi 

(Jy km s" 1 ) 
(7) 

iV H i 

(10 18 cm -2 ) 
(8) 

M H i 

(10 10 Mq) 

(9) 

PG 0007+106 

-211 ± 11 

1054 ± 33 

917 + 22 

5.74 

7.41 

3.81 + 0.14 

13.86 + 0.50 

12.53 + 0.45 

PG 0050+123 

184 ± 21 

1367 ± 64 

916 ± 42 

8.32 

6.80 

4.53 + 0.24 

16.48 + 0.87 

6.40 + 0.34 

PG 0844+249 

120 + 14 

921 ± 43 

474 ± 29 

10.65 

10.06 

4.36 + 0.17 

15.86 + 0.62 

7.61 + 0.30 

PG 1119+120 

-15 + 3 

232 + 8 

222 + 6 

6.23 

7.43 

0.90 + 0.07 

3.28 + 0.24 

1.02 + 0.08 

PG 1211+143 

-235 + 17 

1802 ± 50 

655 ± 34 

13.10 

13.97 

7.89 + 0.21 

28.73 ± 0.78 

22.61 + 0.67 

PG 1351+640 

-332 ± 29 

1202 ± 88 

807 ± 59 

3.56 

4.72 

1.34 + 0.14 

4.89 + 0.52 

4.48 + 0.47 

PG 1411+442 

-205 ± 14 

208 ± 41 

169 ± 27 

3.22 

3.19 

0.42 + 0.08 

1.52 + 0.29 

1.45 + 0.28 

PG 1448+273 

-68 + 12 

168 ± 37 

91 + 25 

3.07 

4.83 

0.26 + 0.04 

0.93 + 0.16 

0.48 + 0.08 

PG 1501+106 

61 + 5 

135 ± 16 

122 + 11 

3.22 

4.53 

0.22 + 0.04 

0.81 + 0.16 

0.14 + 0.03 

F15250+3608 

-89 + 10 

1151 + 30 

1122 + 20 

5.44 

3.67 

2.38 + 0.26 

8.66 + 0.96 

3.23 + 0.36 

PG 1617+175 (U) 

-88 + 13 

1315 + 38 

1208 ± 25 

9.23 

5.81 

6.10 + 0.32 

22.19 + 1.16 

33.32 + 1.75 

PG 2130+099 

67 + 20 

909 ± 61 

366 ± 40 

8.56 

7.87 

2.62 + 0.17 

9.54 + 0.63 

4.37 + 0.29 


Note. — Col.(l): Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic velocity presented 
in Table 1, and its associated measurement error. Col. (3): Velocity width of the emission measured at the 20% level and its error. 
Col. (4): Velocity width of the emission measured at the 50% level and its error. Col. (5): The maximum intensity of the spectrum. 
Col. (6): Signal-to- noise ratio given by the peak intensity relative to the channel rms listed in Table 1. Col. (7): Velocity-integrated 
HI flux in emission and its error. Col. (8): Derived HI column density calculated using Equation 1. Col. (9): Derived HI gas mass 
from emission and its error, calculated using Equation 2. 


TABLE 5 

Measurements of HI in Absorption from Complex Spectral Modeling 


Source 

a) 

^0 

(km s x ) 
(2) 

FWHM 
(km s -1 ) 

(3) 

“Peak” 

(mJy) 

(4) 

S/N 

(5) 

Nhi 

(10 22 cm - 2 ) 
(6) 

Tmax 

(7) 

jabs 

(Jy km s -1 ) 
(8) 

PG 0007+106 

-270 ± 18 

668 ± 43 

-3.26 + 0.18 

3.28 

4.32 

0.033 

-0.38 

PG 0050+124 

-32 + 22 

894 ± 58 

-4.63 + 0.23 

3.41 

112.44 

0.722 

-0.76 

PG 0844+249 

158 ± 22 

700 ± 51 

-4.52 + 0.29 

2.42 



-0.58 

PG 1119+120 

-82 + 2 

25 + 4 

-8.49 ± 1.25 

4.73 



-0.04 


17 + 27 

244 ± 56 

-2.00 + 0.40 

1.11 

29.66 

0.693 

-0.09 

PG 1211+143 

-270 + 4 

466 ± 14 

-37.80 + 1.88 

20.43 

5.82 

0.065 

-3.11 


-61 ± 12 

870 ± 12 

-36.10 + 1.77 

19.52 

10.38 

0.062 

-5.56 

PG 1351+640 

-68 + 18 

952 ± 45 

-3.77 + 0.15 

3.13 

21.32 

0.118 

-0.63 

PG 1411+442 

-288+ 11 

735 ± 27 

-6.19 + 0.20 

5.21 



-0.79 

PG 1448+273 

-106 + 29 

723 ± 69 

-1.98 + 0.16 

1.54 

128.80 

1.079 

-0.26 

PG 1501+106 

-140 + 27 

550 ± 64 

-2.15 + 0.22 

1.47 

>155.41 

1.966 

-0.23 

F15250+3608 

-71 + 21 

858 ± 52 

-5.76 + 0.29 

2.73 

60.56 

0.386 

-0.92 

PG 1617+175 (U) 

-95 ± 30 

870 ± 71 

-4.60 + 0.32 

2.08 



-0.67 

PG 2130+099 

187 ± 35 

325 ± 87 

-2.52 + 0.55 

0.83 

26.35 

0.446 

-0.15 


-247 ± 14 

172 ± 33 

-4.68 ± 0.74 

1.55 

31.27 

1.104 

-0.15 


Note. — Col.(l): Galaxy name. Col. (2): Heliocentric velocity measured from HI relative to the systematic 
velocity presented in Table 1, as modeled using a Gaussian fit, and its associated measurement error. Col. (3): Full 
width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. Col. (4): 
Peak absorption flux and its error as measured using a Gaussian model. Col. (5): Signal-to-noise ratio measured as a 
ratio between the peak flux and channel rms listed in Table 1. Col. (6): Integrated HI column density based on the 
absorption model, assuming a spin temperature of 1000 K, and using Equation 4 based on the optical depth profile 
derived using Equation 3. For PG 1501+106, the 1.4 GHz core flux density is only an upper limit, so its derived 
Ah i value is a lower limit. Col. (7): Peak optical depth based on Equation 3. Col. (8): Velocity-integrated HI flux in 
absorption measured from the best-fit Gaussian component models. 
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TABLE 6 

Measurements of HI Outflow 


Source 

(i) 

Vq 

(km s _1 ) 
(2) 

+ 0, HI 

(km s -1 ) 
(3) 

Nm 

(10 21 cm ~ 2 ) 

(4) 

FWZI 
(km s -1 ) 
(5) 

Vmax 
(km s _1 ) 
(6) 

V 98% M (wide) M (jet) ) 

(km s" 1 ) (M e yr- 1 ) (M 0 yr” 1 ) 

(7) (8) (9) (10) 

ULIRGs 

F04103-2838 (M) 

-409 

-538 ± 36 

81 

1267 + 212 

-1171 ± 214 

-920 ± 70 

1089 

4.4 

41.6 

F05189-2524 

-374 

-374 ± 9 

10 

336 ± 67 

-542 ± 68 

-475 ± 21 

97 

0.4 

40.3 


-79 

-79 ± 13 

12 

459 ± 90 

-309 ± 91 

-218 + 29 

24 

0.1 

38.3 








121+ 

0.5+ 


Mrk 231 

-280 

-280 ± 9 

2 

297 ± 73 

-428 ± 74 

-361 ± 20 

13 

0.1 

39.1 

F15130-1958 (M) 

-847 

-847 ± 10 

39 

266 ± 76 

-980 ± 77 

-923 ± 22 

818 

3.3 

41.9 


-451 

-451 ± 15 

61 

494 + 114 

-697 + 115 

-594 ± 34 

687 

2.7 

41.3 








1505+ 

6.0+ 


F15250+3608 

-71 

18 ± 23 

606 

2744 ± 172 

-1354+ 173 

-53 ± 28 




Median 


39 

459 ± 90 

-697 + 91 

-475 ± 28 

605 

1.6 

40.8 

FIR-strong Quasars 

PG 0050+124 

-32 

-216 + 30 

1124 

2815 ± 191 

-1623 + 192 

-977 ± 58 

6072 

24.3 

41.6 

PG 1119+120 

-82 

-67 + 4 


85 ± 17 

-110 + 17 

-89 + 5 





17 

32 ± 27 

297 

721 ± 182 

-329 ± 184 

-176 + 55 




PG 1211+143 

-270 

-270 + 17 

58 

1680 ± 54 

-1109 + 54 

-666 + 21 

392 

1.6 

40.6 


-61 

-61 ± 21 

104 

3118 ± 55 

-1620 + 56 

-801 ± 23 

158 

0.6 

38.9 








550t 

2.2+ 


PG 1244+026 

-381 

-381 ± 18 

152 

696 ± 62 

-729 ± 63 

-573 ± 25 

1450 

5.8 

41.4 

PG 1351+640 

-68 

264 ± 34 

213 

2946 ± 147 

-1210 ± 148 

-547 ± 51 




PG 1440+356 

325 

325 + 4 


1104 + 45 

-227 ± 45 

41 + 8 




PG 1501+106 

-140 

-201 ± 27 

> 1554 

1637 ± 199 

-1020 + 201 

-669 + 61 

> 7813 

> 31.3 

> 41.6 

PG 2130+099 

-247 

-314 + 24 

313 

547+ 111 

-587 + 112 

-460 ± 37 

2453 

9.8 

41.5 

Median 


213 

1370 ± 87 

-1020 + 88 

-560 + 31 

1450 

5.8 

41.5 

FIR-weak Quasars 

PG 0007+106 

-270 

-59 ± 21 

43 

2053 ± 133 

-1085 + 134 

-627 ± 42 

64 

0.3 

38.5 

PG 0844+249 

158 

38 + 26 


2198 ± 174 

-1061 ± 175 

-558 + 51 




PG 1411+442 

-288 

-83 ± 18 


2362 ± 91 

-1265 + 92 

-709 ± 29 




PG 1448+273 

-106 

-38 + 31 

1288 

2129 + 215 

-1102 + 217 

-653 ± 67 




PG 1617+175 (U) 

-95 

-7 + 33 


2737 ± 234 

-1376 + 236 

-748 ± 69 




Median 


666 

2198 ± 174 

-1102 ± 175 

-653 + 51 

64 

0.3 

38.5 


Note. — * The sum of the mass outflow rates from the multiple Gaussian components in the galaxy. Col.(l): Galaxy name. Col. (2): Heliocentric 

velocity of the absorption component as tabulated in Tables 3 and 5. Col. (3): With the exception of PG 1211+143 and PG 1440+356, this is the outflow 
velocity measured relative to the HI systemic velocity measured from HI emission using values tabulated in column (2) of Table 2 or 4. V©, hi = V© if no 
emission is detected. For PG 1211+143 and PG 1440+356, this velocity is relative to the systemic velocity measured from optical emission lines, since the 
detection of HI emission is uncertain (Figures 5 and 6). Col. (4): Integrated HI column density based on the absorption model from Tables 3 and 5. Col. 
(5): Full- width at zero intensity based on the best-fit Gaussian absorption model and its associated error. Col. (6): The maximum outflow velocity where 
Vmax = V©, hi — ^FWZI. Col. (7): The velocity at which 98% of the gas has lower velocity where Vgg% = V©, hi — 2( )• Col. (8): Derived mass 

outflow rate calculated using Equation 5 and assuming the outflow is in a wide-angled cone where the mass is flowing into a solid angle of n from a radius 
of 1 kpc following the same assumptions as Morganti et al. (2005b). Col. (9): Same as column (8) but for a jet whose opening angle is assumed to be 5° 
launched at a distance of 0.5 kpc. Col. (10): Log of the kinetic energy outflow rate as calculated using Equation 6 for the jet scenario. 
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Fig. 1. — Full intensity GBT spectra of 25 galaxies whose velocities did not overlap with RFI signals. Columns from left to right are 
ordered by source type: ULIRGs - FIR-strong quasars - FIR-weak quasars. Each column is organized with increasing infrared luminosity 
going down. The targets with marginal or uncertain detections as noted in Section 3.2 are labeled with M and U, respectively, after 
their names. The x-axis is displayed as velocity relative to the heliocentric radial velocity of the galaxy given in Table 1 and identified 
by the vertical dashed line. The regions with hash marks are frequencies where the data are affected by persistent RFI spikes in some 
observing sessions, and the black regions are those where the data are affected by RFI spikes in all observing sessions. The best Gaussian 
fits to the absorption features are highlighted in blue, where the dotted lines are individual Gaussian components (when there are multiple 
overlapping components) and the solid curves are the sums of these components. The kinematics of these systems change from narrow 
absorption features likely due to outflows supported by the starbursts in ULIRGs, to broad absorption features from jet-driven outflows 
supported by the AGN in FIR-strong quasars, and finally to the emission profiles more characteristic of gas in rotation in the FIR-weak 
quasars. 
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Fig. 2. — The distribution of equivalent widths of the emission and absorption components as listed in Tables 2 and 3. The absorption 
components are shown as having negative EWs. Targets with only upper limits to their continuum flux densities are excluded from this 
plot. Along the merger sequence, the ULIRGs have the weakest features in both emission and absorption, the FIR-strong quasars have the 
strongest absorption features with moderate strength emission features, and the FIR-weak quasars are dominated by the strongest emission 
features. 
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Fig. 3. — The velocity widths of the components measured in emission (top; W50 in Table 2) and absorption (bottom; FWHM in Table 3). 
The colors represent ULIRGs (black), FIR-strong quasars (red), and FIR-weak quasars (blue). The width of the HI emission in the top panel 
may be inversely correlated with the SFR. The galaxies with the weakest star formation (FIR-weak quasars) have the strongest emission 
features, while the galaxies with the strongest star formation (ULIRGs) have the narrowest velocity widths. The velocity widths of the 
FIR-strong quasars fall between those of the ULIRGs and the FIR-weak quasars. There also seems to be a weak correlation between the 
FWHM in absorption and the SFR for ULIRGs, but this result is very uncertain as ULIRGs tend to have multiple absorption components 
for a single SFR. There is no apparent correlation between the SFR and the absorption parameters in the quasars, suggesting that the 
outflow-related absorption components are not the result of star formation. 
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Fig. 4. — The spectra of the 13 detected galaxies that show polarized HI spectra. The organization of this figure is the same as that of 
Figure 1. The colors represent the two different polarizations (red: XX; blue: YY). Polarized spectra are clearly seen in PG 0007+106, PG 
0050+124, PG 1119+120, PG 1211+143, PG 1351+640, PG 1411+442, PG 1440+356, and F15250+3608. Weak polarized features may 
be present in PG 0844+349, PG 1448+273, PG 1501+106, PG 1617+175 (U), and PG 2130+099. The FIR-strong quasars show the most 
distinct polarized features, perhaps reflecting the presence of polarized radio jets in these objects are contributing polarized continuum 
emission against which the absorption features are seen. 
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Fig. 5. — Left: GBT spectra of PG 1211+143 obtained at four epochs separated by at least one month, with the YY spectrum above and 
the XX spectrum below. The shifts in the central velocity of the deep absorption component are denoted by the vertical dotted lines in 
their corresponding colors. The central velocity shifted by ~120 km s“ 1 over the seven-month period covered by our observations. Right: 
The 2011 November spectrum broken down into four 30-minute periods. There is clear variability in the shape of the spectrum, suggesting 
a varying continuum against which the absorption is detected. If the variability is real, it is consistent with the jet-induced outflow picture 
suggested by Morganti (2011). 
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Fig. 6. — Same as Figure 5, but for PG 1440+356. Short timescale variations are seen in this source, similar to those in PG 1211+143. 
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Fig. 7. — Background-subtracted 0.5-10 keV X-ray light curves of PG 1211+143 (left panel) over 6.5 years through four separate 
observations, and of PG 1440+356 (right panel) over one week through three separate pointings, from XMM-Newton. For each epoch, each 
point represents the EPIC-pn count rate in 1800 s intervals starting with the beginning of the observation in the first bin. The dotted lines 
are the median rates for each specific epoch. The periods of high background flares have been removed from the plot and the calculations 
of the median rate. The error bars are 3— cr. The AGN have both high and low flux states and are clearly variable on the 30-minute time 
scale seen in Figures 5 and 6. 
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Fig. 8. — The HI screen and source configuration that explains the variable polarization-dependent spectra observed in PG 1211+143 
and PG 1440+356. In this picture, the HI screen is extended, covering both the core and the radio jets. The jets produce the polarized 
continuum against which the absorption features are detected. As the central source fluctuates, the ratio of central-to-jet emission also 
changes, resulting in the observed spectra. 
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Fig. 9. — Complex decomposition of the polarized HI spectra assuming that the difference in the two polarizations is due to more HI 
absorption in one polarization than the other. For each galaxy, the upper panel is the difference spectrum between the two polarizations, 
i.e., the spectrum in the polarization with the more negative flux density (Poll) minus the spectrum in the polarization with the higher 
flux density (P 0 I 2 ), plotted in purple. The black curve shows the Gaussian model component(s) best fitting to the difference spectrum as 
tabulated in Table 5. The lower panel shows the residuals of the best fit model as applied to the Poll spectrum and the P 0 I 2 spectrum to 
demonstrate the quality of the fit. The colors in all of the bottom panels, the x-axis, and the hashed rectangles are the same as those in 
Figure 4. The polarization dependence of the HI spectra can be explained by broad HI absorption affecting one polarization more than the 
other. 
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Fig. 10. — Same as Figure 2 apart from the use of the complex fitting results when available (Tables 4 and 5). The general trends with 
stage along the merger sequence are in agreement with those in Figure 2. 
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Fig. 11. — Same as Figure 3 apart from the use of the complex fitting results in Section 5.2 when available (Tables 4 and 5). The general 
trends are in agreement with those in Figure 3. 
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Fig. 12. — The distribution of FWZI of the absorption components in galaxies with detected outflows with Vq hi less than —50 km s 1 
(see Table 6). There is a tentative progression of velocity widths from ULIRGs to FIR-strong quasars to FIR- weak quasars. 


